Fe implanted ZnO: magnetic precipitates versus dilution 
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Nowadays ferromagnetism is often found in potential diluted magnetic semiconductor systems. 
However, many authors argue that the observed ferromagnetism stems from ferromagnetic precipi- 
tates or spinodal decomposition rather than from carrier mediated magnetic impurities, as required 
for a diluted magnetic semiconductor. In the present paper we answer this question for Fe-implanted 
ZnO single crystals comprehensively. Different implantation fiuences and temperatures and post- 
implantation annealing temperatures have been chosen in order to evaluate the structural and 
magnetic properties over a wide range of parameters. Three different regimes with respect to the Fe 
concentration and the process temperature are found: 1) Disperse Fe^""" and Fe^"*" at low Fe concen- 
trations and low processing temperatures, 2) FeZn2 04 at very high processing temperatures and 3) 
an intermediate regime with a co-existence of metallic Fe (Fe°) and ionic Fe (Fe^+ and Fe^+). Ferro- 
magnetism is only observed in the latter two cases, where inverted ZnFe204 and a-Fe nanocrystals 
are the origin of the observed ferromagnetic behavior, respectively. The ionic Fe in the last case 
could contribute to a carrier mediated coupling. However, their separation is too large to couple 
ferromagnetically due to the lack of p-type carrier. For comparison investigations of Fe-implanted 
epitaxial ZnO thin films are presented. 
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I. INTRODUCTION 



Recently, considerable interest has been paid to "spin- 
tronics" , where the spin degree of freedom of the elec- 
tron carries information in the device. One of the ma- 
terial systems to realize this function are diluted mag- 
netic semiconductors (DMS). In DMS materials, transi- 
tion or rare earth metal ions are substituted onto cation 
sites and are coupled with the free carriers to yield fer- 
romagnetism via indirect interaction. In 2000, Dietl et 
al. []| proposed the mean-field Zener model to under- 
stand the ferromagnetism in DMS materials. It has been 
successfully used to describe the magnetic coupling in 
(Ga,Mn)As and (Zn,Mn)Te materials. This model pre- 
dicts that wide bandgap semiconductors (GaN and ZnO) 
doped with Mn exhibit ordering temperatures above 300 
K, provided that a sufficiently large hole density can be 
achieved (10^° cm""^). Sato et al. used the Korringa- 
Kohn-Rostoker Green function method based on the lo- 
cal density approximation of density functional theory to 
calculate the properties of n-type ZnO doped with the 
M TM ions (V, Cr, Mn, Fe, Co, and Ni) [|. The ferro- 
magnetic state, with a Tp of around 2000 K, is predicted 
to be favourable for V, Cr, Fe, Co, and Ni in ZnO while 
Mn-doped ZnO is antiferromagnetic. These predictions 
have largely boosted intensive experimental activity on 
transition metal doped GaN and ZnO. A large number 
of research groups reported the experimental observation 



of ferromagnetism in TM (from Sc to Ni) doped ZnO 
li[i[il@0iS[3[n| fabricated by various methods 
including ion implantation. In contrast to these re- 
sults, other groups rep orted the observations of anti- 



ferromagnetism |12l | [1 
and paramagnetism 



Wm 



spin-glass behavior [Tsj p^. 



18| in TM-doped ZnO. Re- 
cently it was also found that nanoscale precipitates can 
substantially contribute to the ferromagnetic properties 
ildllilKallllKlildKl^ll^. Onc method to intro- 
duce magnetic dopants into ZnO is ion implantation. It 
has several advantages, namely the reproducibility, the 
precise control of the ion fiuence, the use of an isotopi- 
cally pure beam, and the possibility to overcome the sol- 
ubility limit. The major drawback of ion implantation 
is the generation of structural defects in the host lattice. 
However, in several studies it has been demonstrated that 
ZnO exhibits a high amorphization threshold. Therefore 
ion implantation is widely used to dope ZnO with tran- 
sition metal ions. Ref. [28[ gives a review on transition 
metal ion implantation into ZnO. 

Hydrothermal growth is one of the major methods to 
fabricate high quality ZnO single crystals [l^jl^]- Hy- 
drothermal grown ZnO single crystals have been widely 
used for photodiodes [HI, light emitters [s^l, DMS 
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j23j |33j [3411 , and substrates for homoepitaxial growth of 
ZnO films tH. 

The present paper is dedicated to a comprehensive in- 
vestigation of the structural and magnetic properties of 
Fe implanted ZnO bulk crystals grown by the hydrother- 
mal method. Different implantation parameters, i.e. ion 
fiuence, ion energy, and implantation temperature, were 
varied. Three charge states or occupied sites, i.e. metallic 
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Fe and monodispersed Fe^+ and Fe"^"'" and Fe^+ in Zn- 
ferrites, are identified. Metallic Fe nanocrystals (NCs) 
form after implantation at high fluence and high temper- 
ature. They are the major contribution to the measured 
ferromagnetism. The difference between ZnO single crys- 
tals and epitaxial thin films upon the same implantation, 
and the difference between the high (623 K) and the low 
(253 K) temperature implanted samples subject to the 
same annealing, will be discussed. 

Actually the phase separation, namely MnAs precip- 
itates, in (Ga,Mn)As (the most well understood DMS 
material), has been intensively investigated [s^. Fer- 
romagnetic MnAs precipitates are epitaxially embed- 
ded inside the GaAs matrix, and exhibit interesting 
magneto-transport properties 37 1 111 [3§| [i^l . Sato and 
Katayama-Yoshida [4ll | [4^ calculated the chemical pair 
interaction between two magnetic impurities in DMS ma- 
terials. A strong attractive interaction between magnetic 
impurities has been found, which accelerates the spin- 
odal nano-decomposition under thermal non-equilibrium 
crystal growth conditions. Moreover the attractive inter- 
action in (Ga,Mn)As is one order of magnitude weaker 
than that in (Ga.Mn)N (wide bandgap semiconductor). 
Therefore, phase separation in TM doped wide bandgap 
semiconductors is highly expected. In addition to 
those investigations of MnAs/GaAs hybrids, some inter- 
esting magneto-transport properties are demonstrated, 
e.g. anomalous hall effect, and giant magnetoresistance, 
for magnetic NCs embedded inside ZnS [11] and Ge 
[i^ . respectively. Therefore, Fe NCs embedded inside 
ZnO, which are granular magnetic nano-precipitates in- 
side a semiconductor (granular magnetic semiconductor, 
CMS), could have some potential applications in future 
nano-spintronics [i^fisj. 

The paper is organized as follows. First, all the ex- 
perimental methods employed will be described. Then 
the results will be separated according to the physical 
phenomena as follows: lattice damage and recovering, 
the distribution of implanted Fe, the formation of pre- 
cipitates (metallic Fe, or Zn-ferrites) , the charge state of 
Fe, the ferromagnetic properties, and the Fe implanted 
epitaxial ZnO films. In the discussion part, we sketch a 
phase diagram of Fc in ZnO, and apply a model to explain 
the Fe nanocrystal aggregation. Moreover the reason for 
the absence of ferromagnetism in ionic Fe diluted ZnO is 
discussed. 



II. EXPERIMENTS 

Commercial ZnO bulk crystals were implanted with 
^^Fe ions at temperatures ranging from 253 K to 623 K 
with fluences from 0.1 x 10^^ cm~^ to 8 x 10^^ cm~^. 
The implantation energy was 180 keV, which results in 
a projected range of Rp = 89 ± 29 nm, and a maximum 
atomic concentration from 0.14% to 11% (TRIM code 
[46|). For comparison, epitaxial ZnO thin films grown 
on AI2O3 by pulsed laser deposition were implanted with 



^^Fe at selected implantation parameters (623 K, 4 x 10^^ 
cm~^). Three sample series are investigated and listed 
in Table U 

The lattice damage induced by implantation was evalu- 
ated by Rutherford backscattering/channeling spectrom- 
etry (RBS/C). The RBS/C spectra were collected with 
a collimated 1.7 MeV He"*" beam at a backscattering an- 
gle of 170°. The sample was mounted on a three-axis 
goniometer with a precision of 0.01°. The channeling 
spectra were collected by aligning the sample to make 
the impinging He"*" beam parallel with the ZnO<0001> 
axis. Xmin is the channeling minimum yield in RBS/C, 
which is the ratio of the backscattering yield at chan- 
neling condition to that for a random beam incidence 
[IJ. Therefore, Xmm labels the lattice disordering de- 
gree upon implantation. An amorphous sample exhibits 
a Xmin of 100%, while a perfect single crystal corresponds 

to a Xmin 

of 1-2%. 

The Fe distributions were investigated by secondary 
ion mass spectrometry (SIMS), using a Riber MIQ-256 
system with oxygen primary ions of 6 kV and monitoring 
positive secondary ions. The depth scale was calibrated 
by measuring the sputtered crater via profilometry. The 
absolute concentration was determined by calculating the 
sensitivity factors from the low fluence implants. 

Structural analysis was achieved both by synchrotron 
radiation x-ray diffraction (SR-XRD) and conventional 
XRD. SR-XRD was performed at the Rosscndorf beam- 
line (BM20) at the ESRF with an x-ray wavelength of 
0.154 nm. 26-9 scans were used to identify crystalline 
precipitates. 

Conversion electron Mossbauer spectroscopy (GEMS) 
in constant-acceleration mode at room temperature (RT) 
was used to investigate the Fe lattice sites, electronic con- 
figuration and corresponding magnetic hyperfine fields. 
The spectra were evaluated with Lorentzian lines using 
a least squares fit [i^. All isomer shifts are given with 
respect to a-Fc at RT. 

The magnetic properties were measured with a super- 
conducting quantum interference device (SQUID, Quan- 
tum Design MPMS) magnetometer in the temperature 
range of 5-350 K. The samples were measured with the 
field aligned cither along the in- or out-of-planc direction. 
The temperature dependence of the magnetization was 
studied at a constant field and the field dependence at a 
constant temperature. By magnetic measurement, virgin 
ZnO is found to be purely diamagnetic with a suscep- 
tibility of -2.65x10"^ emu/Oc-g. This background was 
subtracted from the magnetic data. To measure the tem- 
perature dependent magnetization after zero field cooling 
and field cooling (ZFC/FC), the sample was cooled in 
zero field from above room temperature to 5 K. Then a 
50 Oe field was applied. The ZFC curve was measured 
with increasing temperature from 5 to 300 (or 350) K, 
after which the FC curve was measured in the same field 
from 300 (or 350) to 5 K with decreasing temperature. 
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FIG. 1: (a) Representative RBS random and channeling spec- 
tra of Fe implanted ZnO with the implantation energy of 180 
keV. The fluence is indicated on the channeling spectra. The 
dashed line separates the damage regions of surface and bulk, 
where the number of displaced atoms is maximum, (b) The 
ion fluence dependence of the maximum relative disorder of 
the Zn lattice (xmin) at different depth (surface and bulk). 



III. RESULTS 

In this section, we present experimental data on struc- 
tural and magnetic properties of '^^Fe implanted ZnO. Of 
interest in this study are the ion-implantation induced 
lattice damage, the distribution of Fc, the formation of 
metallic Fe nanocrystals, the charge state of Fe, the mag- 
netic properties, and the structure and magnetism evo- 
lution upon post annealing. The difference between ZnO 
bulk crystals and epitaxial films upon Fe implantation is 
also compared. 



surface Zn atoms. The implanted Fe ions cannot be de- 
tected for the lowest fluence (0.1x10"'^^ cm~^). However 
they are more pronounced as a hump in the random spec- 
trum for a high fluence of 8x10^^ cm~^. The channel- 
ing spectrum of a virgin sample is provided for compari- 
son. The yield increase in the channeling spectra mainly 
originates from the lattice damage due to implantation. 
However, in the higher fluence case, the Fe ions also sig- 
nificantly increase the RBS yields. Two features are ob- 
served in the RBS/C spectra. One is the bimodal [4§| 
distribution of maximum damage depths, i.e. in the 
bulk and at the surface, separated by the dashed line 
in Figure 1(a). Similar depth pro files have already been 
discussed by Kucheyev et al. [40|. In the bulk damage 
region the nuclear energy-loss profile is maximum, which 
induces a large number of atomic displacements. The 
surface damage peak is often a sink for ion implantation 
induced point defects [49| . 

Another feature is the saturation at larger fiuences. 
Xmin, the ration of the channeling spectrum to the ran- 
dom one, is calculated in both damage regions, as shown 
in Figure 1(b). Above a fluence of 2x10^^ cm~^, both 
damage peaks saturate. This is due to the strong dy- 
namic annealing effect, i.e., migration and interaction 
of defects during ion implantation [4§|. This strong dy- 
namic annealing also makes ZnO an irridation-hard ma- 
terial, i.e., it still partly persists a crystalline state after 
irradiation by Fe ions up to a fluence of 8x10^^ cm~^ 

(Xmm of 68%). 

2. Implantation temperature dependence 

In general increasing the substrate temperature during 
implantation can suppress the lattice damage in semicon- 
ductors. However this is not the case for ZnO. Figure 2(a) 
shows the channeling spectra for Fe implanted ZnO at dif- 
ferent implantation temperatures. Although the surface 
damage peak increases drastically with decreasing im- 
plantation temperature, the bulk damage peak is hardly 
effected by implantation temperature. This can be ob- 
served clearly in Figure 2(b). The point defects induced 
by ion-beam can be significantly suppressed by increasing 
the implantation temperature above 623 K. This tem- 
perature is very critical, and below 623 K, the surface 
damage peak also has no dependence on the substrate 
temperature. This is very important for the electrical 
doping of ZnO by ion implantation, where point defects 
are believe to decrease the conductivity [1^ . 



A. Lattice damage accumulation 

1. Fluence dependence 

Figure 1(a) shows representative RBS/C spectra for 
different Fe fluences implanted at 623 K. The arrow la- 
belled Zn indicates the energy for backscattering from 



3. Recovering by post-annealing 

As shown above, the bulk damage cannot be sup- 
pressed by increasing the implantation temperature. It 
has to be removed by post annealing at higher temper- 
ature. The annealing was performed in high vacuum in 
order to avoid extrinsically induced oxidation of Fe. The 
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FIG. 3: (a) Fe distribution in the samples implanted -with Fe 
at 623 K for different fluences measured by SIMS. The ion 
fluences are indicated in the figure, (b) Distribution of Fe 
in ZnO implanted at different temperatures. The TRIM [i^ 
simulation is presented for a comparison. 



FIG. 2: (a) Representative RBS/C spectra with different im- 
plantation temperature. The fluence is 4xf0^® cm^'^, and 
implantation energy is 180 keV. (b) The calculated Xmin for 
different implantation temperature, Implantation at low tem- 
perature (<473 K) results in more damage at the surface re- 
gion. 



temperature was varied from 823 K to 1073 K. The de- 
tails have been reported in the Ref . [131 ■ Both the surface 
and bulk damage peaks decreased progressively with in- 
creasing the annealing temperature and time. However 
even after annealing at 1073 K for 3.5 hours, there is still 
considerable damage. This is because of the high melting 
point of ZnO (^2250 K). The extended defects can only 
be removed completely by annealing at approximately 
two-thirds of the melting temperature pl\ . Therefore 
a high annealing temperature (1500 K) is necessary to 
completely recover the lattice structure of ZnO. How- 
ever, high vacuum annealing above 1000 K also lead the 
decomposition of ZnO [s^ . 



is smaller than Zn, it is difficult to obtain the depth pro- 
file of the implanted Fe. Therefore, SIMS is employed 
to determine the Fe depth profile (see Figure 3(a)). It 
is observed that the peak concentration of Fe increases 
from 0.1% to 7%, with a projected range of Rp = {80- 
90)± (20-30). This is in a rather good agreement with 
TRIM simulations [i^ . The only discrepancy is the high 
flucncc sample (4x 10^^ cm~^), where TRIM simulations 
predict a peak concentration of 5%. This is due to the 
change in SIMS sensitivities for different materials, which 
implies that the determined concentrations are more ac- 
curate in the low concentration (below 1%) regime. 

In Figure 3(b) the Fe depth profile is compared for 
different temperatures with the same fluence of 0.4 x 10^^ 
cm^'^. The profile does not change significantly due to 
elevating the implantation temperature from 253 K to 
623 K. The slightly higher concentration for implantation 
at 253 K is within the fluence error. 



B. Fe distribution 

RBS/C can give an overview of the lattice damage 
upon Fe implantation. However, since the mass of Fe 



As discussed the in Ref. [l^], Fe diffuses towards the 
surface after high temperature annealing. The same dif- 
fusion of Fe upon annealing was also observed by SIMS 
(not shown). 
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FIG. 4: SR-XRD 20-61 scans of Fe implanted ZnO for diflFerent 
fluences reveal the formation of crystalline Fe nanoparticles. 



Formation of Fe NCs 
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FIG. 5: SR-XRD 26-6 scans of Fe implanted ZnO with the 
same fluence of 4x10^® cm~^ at parts. Only the sample im- 
planted at 623 K shows a-Fe precipitates. 



By employing SR-XRD, we have systematically inves- 
tigated the formation of Fe NCs, and its dependence on 
the fluence and implantation temperature by SR-XRD. 



be crystallographically oriented inside ZnO matrix. This 
highly ordered orientation allows them to be detected 
even by laboratory XRD [2i|[5i|. 



1. Fluence dependence 

Figure 4 shows the SR-XRD pattern (focused on 
Fe(llO) peak) as a function of fluence. At a low flu- 
ence (0.1x10^6 to 0.8x10^6 cm-2), no crystalline Fe NCs 
could be detected, while above a fluence of 2x 10^^ cm^^), 
an Fe(llO) peak appears and increases with fluence. The 
inset shows a wide range scan for the high fluence sam- 
ple (4x1016 cm-2). No other Fe-oxide (FeaOg, Fe304, 
and ZnFe204) particles are detected in the as- implanted 
state. The full width at half maximum (FWHM) of the 
Fe(llO) peak decreases with fluence, indicating a growth 
of the average diameter of these NCs. The crystallite size 
is calculated using the Scherrer formula [ssj . 

(i=0.9A/(/3-cos6') (1) 

where A is the wavelength of the x-ray, 9 the Bragg angle, 
and (3 the FWHM of 29 in radians. The crystallite size 
was estimated using Eq. [T] and is listed in table HIl 

Note, that only one peak of Fe(llO) appears in the in- 
set of Figure 4. This indicates a texture of the Fe NCs. 
However no texture behavior is found even for the high- 
est fluence sample in pole figure measurements on Fe(llO) 
and Fe(200) (not shown). This could be due to the differ- 
ence in the crystalline symmetry of hexagonal ZnO (six 
fold symmetry) and bcc-Fe (four fold symmetry). For a 
bcc-crystal, one cannot find a six-fold symmetry viewed 
from any direction. In contrast, hcp-Co(OOOl) and fcc- 
Ni(lll) NCs, which are six-fold symmetric, are found to 



2. Implantation temperature dependence 

SR-XRD was also performed for the samples with an 
Fe fluence of 4x 10^^ cm~^ implanted at different temper- 
atures from 253 K to 623 K. As shown in Figure 5, for 
implantation temperatures of 473 K and below, no crys- 
talline Fe could be detected. This is also conflrmed by the 
CEMS results (shown later), where the Fe" state appears 
only at an implantation temperature of 623 K. Note the 
asymmetry of the ZnO (0002) diffraction peaks in Figure 
5. Shoulders on the right side (smaller lattice constant) 
are clearly observed. These shoulders decrease with in- 
creasing implantation temperature, and can therefore be 
associated with lattice damage or ZnO substituted with 
Fe. In view of a detailed study of ion implantation into 
GaN where the implantation induces a lattice expansion 
of GaN (a shoulder at left side) [sHilEHl, we rather at- 
tribute the observed right side shoulders to ZnO substi- 
tuted with Fe. In the 623 K implantation, metallic Fe 
NCs start to form, therefore the substitution is reduced. 

3. Growth with post-annealing 

After thermal annealing at 823 K for 15 min. more a-Fe 
NCs of a larger size are formed for the samples implanted 
at 623 K. After 1073 K and 15 min annealing, the a- 
Fe almost disappears and ZnFe204 starts to form. The 
details of the structure and magnetism evolution upon 
thermal annealing can be found in the Ref. [l^. 
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TABLE I: Structural properties of ^'^Fe-implanted ZnO bulk crystals and epitaxial thin films. The implantation energy is 180 
keV. 



Fluence 
(cm-2) 


Implant at ion 
Temperature (K) 


Peak Concentration 
(TRIM simulation) 


Peak Concentration 
(SIMS) 


Bulk 


Surface 


Metainc re 
Formation 


OTxlO"" 


623 


0.14% 


0.1% 


8.1% 


5.9% 


No 


0.4x10^^ 


623 


0.55% 


0.46% 


38% 


14% 


No 


0.8x10^*^ 


623 


1.1% 


0.89% 


43% 


16% 


No 


2x10 


623 


2.7% 


2.6%) 


57%) 


39%) 


a-Fe 


4x10^® 


623 


5.5% 


6.0% 


60% 


32% 


a-Fe 


8x10^® 


623 


11% 




65% 


55% 


a-Fe 


0.4x10^^ 


253 


0.55% 


0.5% 


31% 


16% 


No 


4x10^*^ 


253 


5.5% 


5.5% 


65% 


43% 


No 


4x10^'^ 


298 


5.5% 




65% 


40% 


No 


4x10^*^ 


473 


5.5% 




66% 


42% 


No 


4xl0^^(") 


623 K 


5.5% 




44% 




a and 7-Fe 



"ZnO epitaxial thin films. 



TABLE II: Structural and magnetic properties for Fe-implanted ZnO. The ferromagnetic fraction corresponds to the percentage 
of ferromagnetic Fe (at 5 K) compared with all implanted Fe ions. The crystallite size evaluated by ZFC magnetization is only 
for Q-Fe NCs. 





Sample 




Crystallite size 


Tfl 


Crystallite size 


Saturation 


Ferromagnetic 


Coercivity 


Fluence 


Timp. 




(XRD) 


(ZFC) 


EqlU 


Eq. m 


magnetization (5 K) 


Fe fraction 


at 5 K 


(cm-^) 


(K) 


(K) 


(nm) 


(K) 




(nm) 


(/is/Fe) 


(%) 


(Oe) 


2x10"" 


623 




5.6 


38 


8 


6.6 


0.08 


3.6 


600 


4x10^^ 


623 




7.1 


137 


12 


8.9 


0.24 


11 


360 


8x10^® 


623 




8.9 


212 


14 


11.3 


0.13 


5.9 


360 


4x10^^ 


623 


823 


9.4 


200 


14 


10.2 


0.34 


15 


360 


4x10^*^ 


253 


823 




295 


16 


9.5 


0.52 


24 


370 


4x10^^" 


623 




-/6 (a/7-Fe) 


26 


7 


4.6 


0.55 


25 


220 




623 


823 


8.1/11 (a/7-Fe) 


280 


15 


10.2 


1.3 


59 


220 



"ZnO epitaxial layers implanted with Fe at 180 keV and 623 K. 
''Calculated by comparing the saturation magnetization with the 
value (2.2 Ms /Fo) for bulk Fe. 



D. Charge state of Fe 



Fluence dependence 



GEMS allows one to identify different site occupa- 
tions, charge and magnetic states of ^^Fe. The hyper- 
fine parameters calculated according to the evaluations 
of the spectra are given in Tabic IIIII All isomer shifts 
are given relative to an a-Fe reference foil. In general, 
the implanted Fe occupy three different states: metal- 
lic Fe, Fe^+ and Fe'^^ ions dispersed in the ZnO matrix 
and finally Fe^"*" in Zn-ferrites. The hyperfine interac- 
tion parameters obtained from the best fits are different 
from that of ferromagnetic a or 7-Fe203 or ferrimagnctic 
FC3O4. Hence, the presence of these phases was excluded. 
The dispersed ionic Fc could substitute onto Zn site. 



Figure 6(a) and 6(b) show the comparison of Fe im- 
planted ZnO at 623 K with a fluence of 0.4x10^'^ and 
4x10^® cm^^, respectively. In spectrum (a), the sin- 
glet S and doublet D(I) arc attributed to Fe^+, while 
the doublet D(II) is from Fe^+. In the high fluence sam- 
ple (spectrum (b)), the majority of Fe arc ionic states 
Fc^+ (singlet S) and Fc'^+ (doublet D(I) and D(II)), while 
a considerable fraction of a sextet associated to a-Fe is 
present (sextet M). The formation of a-Fe is in agree- 
ment with SR-XRD observation (Figure 4). At room 
temperature, all Fe^^ and Fe^^ show no ferromagnetic 
interaction. Later on in subsection IIII Fl we show that 
even at 5 K the measured fcrromagnctism can only be 
attributed to a-Fe NCs. 
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TABLE III: Hyperfine parameters obtained from the evaluation of CEMS for samples implanted or annealed at different 
temperatures. The fluence was 0.4x10^'' cm^^ for the first sample, while 4x10^^ cm^^ for all other samples. The codes of SI, 
S2 and S3 notate the samples for post-annealing process. The notations for the fitting lines are given as S (singlet), D (doublet) 
and M (sextet). 



Code 


Sample 
(K) 


(K) 


FR" 
(%) 
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(%) 


D(I) 
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{mm/s) 
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{mm/s) 
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FR" 

(%) 


M (a-Fe) 

IS' 
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0.57 
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0.31 


0.75 


58.6 


0.81 
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SI 


623 




32.8 


0.53 


31.5 


0.78 


1.29 


23.2 


0.96 


0.58 


12.5 


0.06 


30.5 


SI 


623 


823 


42.6 


0.42 
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1.52 
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0.94 
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18.2 
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31.7 


SI 


623 
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0.32 
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46.8 
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0.88 
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Dist.^ 


S3f 
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26.3 


0.45 


0.32 


37.3 


0.91 


0.80 


13.4 





Dist.= 


S3 


623 


823 


13.9 


-0.09 


31.4 


0.45 


0.35 


15.9 


0.92 


0.58 


38.8 


0.02 


Dist." 


"Fraction corresponding to the relative area 


of the subspcctrum. 

















'Isomer shift: mm/s for a-Fe, 0.7-1.2 mm/s for Fe2+, 0.2-0.7 
mm/s for Fe'^+, and -0.1 mm/s for 7-Fe. 
'^Quadrupole splitting. 
''Magnetic hyperfine field. 
■^Hyperfine field distribution. 
■'^ZnO epitaxial thin films. 
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FIG. 6: Room temperature CEMS for ZnO bulk crystals 
implanted with ^^Fe with different fluences, and with post- 
annealing. The notations for the fitting lines are given as S 
(singlet), D (doublet) and M (sextet). The fluence and the 
process temperatures are indicated. 



2. Implantation temperature dependence 

Figure 7(a) and 7(b) shows CEMS for the samples im- 
planted at low temperatures, 473 K and 253 K, respec- 
tively, with a '^^Fe fluence of 4x10^^ cm~^. In these two 
samples, ionic Fe are the dominant charge states: Fe'^"'" 
(S and D(I)), and Fc2+ (D(II)). In contrast to Figure 
6(b), there is no detectable a-Fe in these two samples. 
This is also in agreement with SR-XRD results (Figure 
5), where up to an implantation temperature of 473 K 
no a-Fe is found. 



E. Implantation Energy dependence 

1. Evolution with post-annealing 

The post-annealing was performed on selected sam- 
ples: SI and S2 (Table [111)) . They were implanted with 
the same fluence of 4x10^^ cm'^ at 623 K and 253 K, 
respectively. For sample SI, upon annealing at 823 K 
for 15 min, the intensity of the sextet increases up to 
18.2% while the fraction of Fe^+ (doublet D(I)) decreases, 
suggesting the growth of the a-Fe nanoparticles and the 
recovery of lattice defects (Figure 6(c)). Moreover, the 
value for the magnetic hyperfine field B?j/ increases upon 
annealing and moves closer (from 30.5 T to 31.7 T) to 
the known value for bulk a-Fe (33 T). For sample S2, af- 
ter annealing at 823 K for 15 min, the relative fraction of 
metallic a-Fe increases up to 28.8%. The hyperfine field 
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FIG. 7: Room temperature GEMS for ZnO implanted with 
^'^Fe at 253 K and subsequent annealed at 823 K for 15 min. 
The notations for the fitting lines are given as S (singlet), D 
(doublet) and M (sextet). On the right side of the spectra, 
the probability distribution P for the magnetic hyperfine field 
(Bhf, solid lines)are given. 



is distributed with maxima at 18 T, 27 T and mostly at 
32.5 T (a-Fc). Comparing with the annealing of sample 
SI, a larger fraction of a-Fe is formed in the annealed 
sample S2. This is consistent with a larger magnetiza- 
tion measured by SQUID (shown later). In addition, a 
small fraction of singlet (S) presents, which is attributed 
to 7-Fe according to the isomer shift. 

Higher temperature (1073 K) annealing was performed 
on sample SI and has been reported in Ref [27j. After 
annealing at 1073 K for 3.5 hours, Fe"^+ is the only charge 
state, and Zn-ferrites (ZnFe204) are formed and are crys- 
tallographically oriented inside ZnO matrix. 



F. Magnetic properties Fe implanted ZnO 

In the previous sections, we have reported a thorough 
investigation on the structural properties, and the charge 
states of Fe. The main conclusion can be summarized as 
follows (i) upon implantation at a temperature of 623 
K, a small part (around 12%) of the implanted Fe ions 
forms as crystalline Fe already in the as-implanted state, 
while the major part of the implanted Fe is in ionic states 
(Fe^"*", and Fe'^^); (ii) implantation at a low temperature 
(253 K) suppresses the metallic Fe formation, and the 
implanted Fe ions are in ionic states, but they are not 
magnetically coupled at room temperature; (iii) post- 
annealing at 823 K largely enhances the Fe NC formation 



in all implanted samples for both implantation temper- 
atures (253 K and 623 K). Since CEMS was performed 
at room temperature only, the magnetic properties of 
metallic and ionic Fe at low temperature could not be 
determined. Here we present the results from SQUID 
magnetometry measured from 5 K to 350 K. We will 
show that the metallic Fe NCs are superparamagnetic, 
aird they are the predomiirant contribution to the mea- 
sured ferromagnetic response even at 5 K. In contrast the 
ionic Fe is not ferromagnetically coupled even at 5 K. 



1. Superparamagnetism of Fe NCs 

For magnetic nanoparticles, the formation of domain 
walls is energetically unfavorable and below a certain size 
(typically in the range of 15 to 30 nm depending on the 
material), the particle stays in a single-domain config- 
uration. The magnetism of a single nanoparticle in a 
solid matrix can be described by the Ncel process [57| . 
If the particle size is sufficiently small, above a partic- 
ular temperature (so-called blocking temperature, T^) 
thermal fluctuations dominate and no preferred magne- 
tization direction can be defined. Such a system of su- 
perparamagnetic particles does not exhibit a hysteresis 
loop above Tb; therefore the coercivity (He) and the re- 
manence (M/j) are both zero. Phenomenologically there 
are two characteristic features in the temperature de- 
pendent magnetization of a nanoparticle system. One 
is the irreversibility of the magnetization in a small ap- 
plied field {e.g. 50 Oe) after zero field cooling and field 
coofing (ZFC/FC) [53|. The other is the drastic drop 
of the coercivity and of the remanence at a temperature 
close to or above Tb [sst . 

For a dc magnetization measurement in a small mag- 
netic field by SQUID, the blocking temperatue T^ is 
given by 



Tp. 



KeffV 

30fcs 



(2) 



where Keff{V) is the anisotropy energy density, V the 
particle volume, and ks the Boltzmann constant [57| . 
With this equation, one cair estimate the particle size 
[5^. However, in any fine particle system, there is a 
distribution of particle sizes, which is usually assumed as 
a log-normal distribution D{V). 



D{V) 



A 



2ti(jlV 



exp[ 



ln[V/V„. 



(3) 



where Vmean is the most probable value, and is the 
standard deviation. 

Such a volume distribution results in a distribution of 
blocking temperatures Tb(V). The ZFC magnetization 
can be calculated as follows [53] ■ 
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FIG. 8: ZFC magnetization vs temperature for the as- 
implanted sample with the fluence of 4x10^® cm~^. Solid 
line is the best fitting using the Eq. |4] Inset shows the size 
distribution of the Fe NCs deduced from the analysis of ZFC 
magnetization. 



where Mg is the spontaneous magnetization of the 
particle, DiV) the volume distribution, VumitiT) (= 
SOksT/Keff) the maximum volume in the superpara- 
magnetic state, N the normalizing factor, and fc^ the 
Boltzmann constant. is assumed to be a constant in- 
dependent of temperature [5^ (soj . The first integral rep- 
resents the contribution of the superparamagnetic parti- 
cles, while the second corresponds to the blocked ones. 
A more precise determination of the size should be per- 
formed by fitting the ZFC curve with the equation U] 
Figure 8 shows the fitting on the ZFC magnetization 
curve, and the corresponding size distribution. In the 
fit. Kpf f{V)^5xlO^ Jm~^ is treated as a constant. 



2. Fluence dependence 

Figure 9(a) shows the ZFC/FC magnetization curves 
in a 50 Oc field for different fluences of Fe implanted ZnO. 
The FC curves for low fluences of O.IxlO^^ (not shown 
for clarity) and 0.8x10^^ cm~^ completely overlap with 
the corresponding ZFC curves at values close to zero. No 
superparamagnetic particles are present in the two sam- 
ples. For larger fluences (above 2x10^^ cm~^), a distinct 
difference in ZFC/FC curves was observed. ZFC curves 
show a gradual increase (deblocking) at low tempera- 
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FIG. 9: (a) Magnetization curves with an applied field of 50 
Oe after ZFC/FC for the Fe implanted ZnO. With increasing 
fluence, the Fe NCs are growing in size, resulting in a higher 
blocking temperature; (b) ZFC/FC curves with different ap- 
plied fields. The blocking temperature decreases progressively 
with increasing field. 



tures, and reach a broad peak with a maximum, while 
FC curves continue to increase with decreasing tempera- 
ture. The broad peak in the ZFC curves is due to the size 
distribution of Fe NCs. In this paper, the temperature at 
the maximum of the ZFC curve is taken as the average 
blocking temperature (later referred as T^). At a much 
higher temperature than T b , FC curves still depart from 
corresponding ZFC curves, which distinguish the Fe par- 
ticle system from a conventional spin-glass system where 
the FC curve merges together with ZFC curve just at 
Tb and shows a plateau below Ts [eO]. The ZFC/FC 
curves are general characteristics of magnetic nanoparti- 
cle systems with a broad size distribution [6l| . T^ in- 
creases with the fluence, i.e. the size of nanoparticles. 
Table [II] lists the average size of Fc NCs calculated by 
Eq. [5] and by XRD data (Eq. [T|), and simulated by Eq. 
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FIG. 10: (a) Hysteresis loops measured at 5 K for Fe im- 
planted ZnO with different fluence. The inset shows the room 
temperature M-H loop for the sample implanted with a flu- 
ence of 4x10^*^ cm^^. (b) Hysteresis loops measured at 5 K 
for Fe implanted ZnO at different implantation temperatures 
with the same fluence of 4x10^^ cm~^. 



m Although the trend is the same for all calculations, 
the values from Eql2] are larger than that from Eq. [4] 
and from XRD data. Given the large size distribution in 
the present magnetic nanoparticle system, is overes- 
timated by ta king the temperature at the maximum of 
the ZFC curve [5^. This explains why Eq. [2]gives a large 
average particle size. For the fitting according to Eq. IH 
one has to note that Mg and K^f f {V) are assumed to be 
temperature independent, and the interaction between 
the NCs is ignored. These effect contribute to the error 
bars [b^]. Nevertheless, from both techniques (XRD, and 
ZFC magnetization), wc have determined the size of Fe 
NCs and its distribution. 

Figure 9(b) shows ZFC/FC curves with an applied field 
ranging from 50 Oe to 500 Oe for the sample with fluence 
of 4x10"^^ cm~^. Tb in the ZFC curves clearly shifts to 
lower temperatures with increasing field. This behavior 
is expected for magnetic nanoparticles since the reduced 
energy barrier caused by the external magnetic field al- 
lows the reorientation of the superparamagnetic moments 
by thermal fluctuations at lower temperatures [1^ [11] . 



Figure 10(a) shows the magnetization versus field re- 
versal (M-H) of samples implanted with large fluences. 
At 5 K, hysteretic behaviors were observed for all three 
samples. The saturation moment is increased with in- 
creasing fluence, however the coercivity is decreased from 
600 Oe for a fluence of 2x10^^ cm^^ 339 Qe for larger 
fluences (see table I). This can be explained by the en- 
hanced coercivity effect for the interfacial spins, which in- 
creases with decreasing the size of nanoparticles . The 
inset shows the M-H curve at 300 K for the sample im- 
planted with the fluence of 4x 10^^ cm~^. As expected for 
a magnetic nanoparticle system, above the blocking tem- 
perature, both remanence and coercivity drop to zero. 



3. Implantation temperature dependence 

Figure 10(b) shows the magnetization versus field re- 
versal of samples implanted with Fe (4x10^^ cm~^) at 
different implantation temperatures. Only the sample 
implanted at 623 K shows a hysteretic behavior due to 
the presence of Fe NCs, while the other samples im- 
planted at 473 K or below show no ferromagnetic re- 
sponse down to 5 K. This is in full agreement with SR- 
XRD and CEMS measurements. 



4- Post annealing effect 

The magnetic properties of the samples implanted at 
623 K have been reported in a previous paper [l^l • The 
main conclusions are the following: the annealing at 823 
K results in the growth of a-Fe NCs. During annealing 
at 1073 K the majority of the metallic Fe is oxidized; af- 
ter a long term annealing at 1073 K, crystallographically 
oriented ZnFe204 NCs form. Here we mainly present 
the annealing at 823 K for the samples implanted at 253 
K. Due to the different initial state from the 623 K im- 
planted samples, the same annealing temperature leads 
to different results. 

Figure 11 shows the magnetic properties of the samples 
implanted at 253 K with subsequent post annealing. In 
the as-implanted state, there is no ferromagnctism down 
to 5 K, while after 823 K annealing, magnetization of 
0.52 TOs/Fe was observed. The ZFC/FC magnetization 
curves show the characteristics of magnetic nanoparticle 
system. According to SR-XRD and CEMS results, we 
attribute this to Fe NCs. The ZFC curve is very broad 
and Tb is above room temperature. The M-H curves at 
300 K for both cases are still open, although with much 
smaller coercivity and remanence compared with 5 K. 
However the magnetic properties are quite different from 
the 623 K implanted sample after post annealing at 823 
K, where the Tb is well below 300 K, and at 300 K there 
is neither coercivity nor remanence. We will discuss this 
difference in section HVl 
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FIG. 11: Hysteresis loops measured at 5 K for ZnO implanted 
with 180 keV Fe at 253 K (as-implanted and post annealed). 
The fluence is 4x10^*" cm~^. For comparison, the sample 
implanted with the same energy and fluence, but at a high 
implantation temperature (623 K), after the same annealing 
process is shown as the solid line, (b) ZFC/FC magnetization 
of the sample after 823 K annealing. The arrows indicate the 
blocking temperatures. 



FIG. 12: Hysteresis loops measured at 5 K for Fe implanted 
ZnO at 180 keV and 623 K up to a fluence of 4x10^'^ 
cm~^. The field is changed from parallel to perpendicular 
with respect to the sample surface, revealing the magnetic 
anisotropy. The intersections between the easy axis M-H 
curve and solid and dashed direct lines indicate the lowest 
and highest Ha- 



difficult to be deduced since it is not easy to measure 
a real hard axis loop with SQUID magnetometry with- 
out a precise control of the sample alignment. Moreover 
the size distribution of Fe NCs could result in a distri- 
bution of Ha- Therefore, we deduce a lower and upper 
limit of Ha according to the shape of the hard axis loop. 
Using this approach, the anisotropy energy is estimated 
to be in the range of (1.8-3.3)x 10^ Jm~^. It is larger 
than the magnetocrystalline anisotropy, and around one 
order of magnitude larger than the uniaxial anisotropy 
observed in Fe thin films [g^l and micro-scale Fe nano- 
magnets If the Fe NCs are assumed to be sphere- 

like, their magnetism should be isotropic, unless they are 
textured. However as found by XRD, these Fe NCs are 
not textured. This magnetic anisotropy could be due to 
the shape effect of Fe NCs, i-e. they are not sphere-like, 
or magnetostriction. There is, however, no evidence for 
any of these two possibilities. 



5. Magnetic anisotropy of Fe NCs 

M-H loops were also measured for selective samples 
which have been implanted with a fluence of 4x10^^ 
cm~^ with the field applied perpendicular to the sample 
surface. Figure 12 shows the comparison of the magneti- 
zation between the in-plane and out-of-plane direction at 
5 K. The in-plane (parallel to the ZnO surface) is the easy 
axis, while the out-of-plane (perpendicular) is the hard 
axis. At 5 K, the coercivity of the easy axis is around 360 
Oe, and the ratio of M^/Mg (remanence and saturation 
moment) is around 58%. The anisotropy energy, K, can 
be calculated according to the equation of K=MaHs/2, 
where Mg is the saturation moment of 47rAfs=22000 G, 
Ha is the effective anisotropy field. Indeed Ha is rather 



6- Memory effect of Fe NCs 

Below the blocking temperature, a magnetic nanopar- 
ticle system has a rich and unusual behavior. For in- 
stance a slow relaxation and a history-dependent mag- 
netic memory are found in the dc magnetization as a 
function of temperature jBC" 



61| [66| [67[ [68| . In our sys- 



tem, Fe nanoparticles embedded inside ZnO crystals, the 
temperature dependent memory effect was also observed 
(Figure 13) using a cooling and heating protocol sug- 
gested by Sun et al. [6^ . At 300 K a magnetic field of 50 
Oe was applied and the sample was cooled down to 5 K at 
a constant cooling rate of 3 K/min. Then the sample was 
heated continuously at the same rate and the magnetiza- 
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FIG. 13: Temperature dependent memory effect in the dc 
magnetization. Tfie reference curve is measured on heating 
at a constant rage of 3 K/min after PC in 50 Oe. The solid 
squires are measured during cooUng in 50 Oe at the same rate 
but with a stop of 2 hours at 50 K. The field is cut off during 
stop. The open circles are measured with continuous heating 
at the same rate after the previous cooling protocol. Inset 
shows the reciprocal magnetization versus temperature. 



tion was recorded. The obtained M(T) curve is referred 
as the reference curve (solid line in Figure 13). There- 
after, we cooled the sample at the same rate and recorded 
the magnetization with cooling, but temporarily stopped 
at T = 50 K for a waiting time of 2 hours. During waiting 
time, the field was switched off. After the stop, the 50 Oe 
field was reapplied and cooling and measuring were re- 
sumed. The temporary stop resulted in a steplike M(T) 
curve (solid squares in Figure 13). After reaching the low- 
est temperature 5 K, the sample was heated back with the 
rate of 3K/min in the same field, and the magnetization 
was recorded again. The M(T) curve during this heating 
also has a steplike behavior at the stop temperature of 50 
K, then recovers the previous M(T) curve measured dur- 
ing cooling. The system remembers its thermal history. 
Two explanations have been suggested for such a mem- 
ory effect The first one is a broad distribution of 
blocking temperatures originating from the distribution 
of the anisotropy energy barriers. Another explanation 
is the strong dipolar interaction between nanoparticles, 
which frustrates the nanomagnetic moments, and slows 
down their relaxation. Our observations rather support 
the first model. First of all, the memory effect is also ob- 
served for two other samples (2x10^^, and 8x10^^ cm~^) 
(not shown). Therefore the effect is independent of ion 
fluence, i.e. particle density. Second, the inset of Fig- 
ure 13 shows the reciprocal FC magnetization at 500 Oe 
versus temperature. The perfect linearity of the curve 
for T>50 K strongly suggests that the dynamics of the 
nanoparticles above blocking temperature can be well de- 
scribed by superparamagnetism. Therefore, the magnetic 
properties of the sample depend only on the individual 
particle behavior. Third, the size of Fe nanoparticles is 
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FIG. 14: XRD 29-9 scans of Fe implanted epitaxial ZnO lay- 
ers: as-implanted and post annealed at 823 K. The bulk crys- 
tal ZnO implanted at the same condition is shown for com- 
parison. 



widely dispersed according to the analysis on the ZFC 
magnetization curve as shown in Figure 8. Therefore, we 
would attribute the memory effect to the broad distribu- 
tion of particle size, i.e. of anisotropy energy barriers. 



G. Fe implanted epitaxial ZnO layers 

The epitaxial ZnO layers used in this study were grown 
by pulsed layer deposition on Al2O3(0001). These thin 
films are n-type conducting with a carrier concentration 
of lO^'^-lO^^ cm^ at room temperature. Details about the 

sample preparation can be found in Refs. [Hlziinl. 

^^Fe ions were implanted at an energy of 180 keV at 
623 K. Then the samples were subjected to the same 
thermal annealing, and structural as well as magnetic 
characterization like the bulk crystals. 



1. Formation of Fe NCs 

Figure 14 shows the XRD 29-0 scans of Fe implanted 
ZnO epitaxial layers (as-implanted and post annealed at 
823 K). In the left panel, one can see the nice epitaxy 
of ZnO on AI2O3 with the out-of-plane relationship of 
ZnO(0001)||Al2O3(0001). The right panel is a zoom on 
the Fe-rclated peak region. A single crystal sample im- 
planted with the same fluence and at the same temper- 
ature is shown for comparison. Obviously, the epitaxial 
ZnO behaves differently from the bulk crystals upon Fe 
implantation. In the epitaxial-layer, 7-Fe is the predom- 
inant phase, while it is a-Fe in the single crystal. This 
difference will be discussed in section IIVI Upon thermal 
annealing at 823 K, the epitaxial ZnO behaves similar 
to the low temperature implanted bulk crystals. Both 
metallic Fe phases (a and 7) are growing. 
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FIG. 15: Room temperature GEMS for ZnO epitaxial thin 
films as-implanted with ^^Fe and post-annealed at 823 K for 
15 min. The notations for the fitting lines are given as S 
(singlet), D (doublet) and M (sextet). On the right side of 
the spectra, the probability distribution P for the magnetic 
hyperfine field (Bhf, solid lines) are given. 



2. Charge state of Fe 

The charge and chemical states of Fe deduced from 
CEMS are shown in Figure 15. The hyperfine interaction 
parameters arc given in Table IIIII In the as-implanted 
sample, ionic Fe is the predominant phase, while also a- 
and 7-Fe are present (sextet M and singlet S, respec- 
tively). After annealing at 823 K, the fraction of fer- 
romagnetic a-Fe is drastically increased from 13.4% to 
38.8%. In the bulk crystal implanted at the same con- 
dition, there is no 7-Fc neither in as-implanted nor in 
annealed samples. Also the fraction of a-Fe after anneal- 
ing (18.2%) is much lower than that in epitaxial films 
(38.8%). 
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FIG. 16: (a) ZFG/FC magnetization curves at an applied 
field of 50 Oe for the as-implanted and 823 K annealed ZnO 
films; (b) M-H loops at 5 K and 300 K. The bulk crystal ZnO 
implanted at the same condition is shown as solid lines for 
comparison. 



magnetization is increased from O.SS/is/Fe to 1.3/iB/Fe 
at 5 K and from 0.24^B/Fe to l.l^s/Fe at 300 K, respec- 
tively, with annealing. For both samples, the M-H loops 
show no hysteresis at 300 K without coercivity and re- 
manence. Obviously, the annealing behavior is different 
from the single crystal implanted at same temperature of 
623 K, but similar to the single crystal implanted at 253 
K. We will discuss this point in section HVl 



The ZFC /FC magnetization measurements for the as- 
implanted and 823 K annealed samples (Figure 16(a)) 
shows the typical behavior of a magnetic nanoparticle 
system. However T^ in the ZFC curves increases from 
26 K to around 300 K with post annealing. The absolute 
magnetization value per Fe in the ZFC /FC curves for the 
annealed sample is lower than that of the as-implanted 
sample. This is due to the fact that there are more bigger 
Fe NCs after annealing, and the bigger NCs are more dif- 
ficult to be aligned at such a small field of 50 Oe. Figure 
16(b) shows the M-H curves. The coercivity is not sig- 
nificantly changed with annealing, while the saturation 



IV. DISCUSSIONS 

A. Phase diagram of Fe in ZnO 

In section Hill we present the structure and magnetic 
properties of Fe implanted ZnO. The implantation pa- 
rameters, i.e. fluence, energy, temperature, were varied. 
In general, metallic Fe NCs have been formed already in 
the as-implanted state when the implantation tempera- 
ture is above 623 K and the fluence is above 2 x 10^^ 
cm~^. By summarizing all results, a phase diagram of 
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Figure 17 



1200 



1000 
ro 800 

Q. 



_Fe in ZnFe 

2 4 



a-Fe and 

dispersed Fe^* and Fe^* 



Dispersed Fe^*, Fe^* 



600 



w 

^ 400 
O 
g 

°- 200 

1 2 3 4 5 6 

Fe concentration (%) 

FIG. 17; The phase diagram of Fe in ZnO bulk crystals de- 
rived from the data presented in this work. The process tem- 
perature refers to the implantation or annealing temperature. 



Fe in ZnO can be sketched, as shown in Figure 17. Note 
that the materials studied in this research are ZnO bulk 
crystals grown by hydro-thermal method. They are semi- 
insulating in the as-purchased state with n-type carrier 
concentration of 10^^-10^^ cm^'^. The phase diagram will 
likely be different for cpitaxial-ZnO and for p-typc ZnO. 



B. Phase separation depends on the forms of ZnO 

In section Hill we have shown the structural and mag- 
netic properties of ZnO bulk crystals and epitaxial thin 
films implanted at the same temperature and with the 
same Fe fluence. They are obviously different from each 
other (see table . In the bulk crystals, only 11% (in- 
creased to 15% after 823 K post annealing) of the im- 
planted Fe is ferromagnetic and mainly a-Fe, while 25% 
in the epitaxial ZnO films (increased to 59% after 823 
K post annealing). Recently Dictl proposed the self- 
organized growth driven by the charge state of the mag- 
netic impurities [4^[7l|. The energy levels derived from 
the open d shells of transition metals reside usually in the 
bandgap of the host semiconductor. The mid-gap lev- 
els of magnetic impurities trap carriers origination from 
residual impurities or defects. This trapping alters the 
charge state of the magnetic ions and hence affects their 
mutual Coulomb interactions. Therefore, different car- 
riers (electrons or holes, with different concentrations) 
could lead to different interactions {e.g. repulsions and 
attractions) between the implanted transition metal ions, 
and finally result in a different phase separation. Both 
ZnO materials (bulk crystals and epitaxial thin films) 
used in this study are n-type semiconductors. The carrier 
concentration is around 10^^-10^'^ cm~^ for bulk crystals 
and 10^^-10^^ cm^^ for epitaxial thin films [6l][7^ at room 



temperature. Therefore, we can explain the different be- 
havior in ZnO bulk crystals and epitaxial layers upon Fe 
implantation in the above-mentioned model. A higher 
concentration of free electrons leads to more agglomera- 
tions of Fe. Moreover the nanocrystal aggregation could 
be largely reduced or avoided by the realization of p-type 
doping in ZnO. 

However, one has to note that there are a lot of defects, 
such as dislocations, and stacking faults, in the epitax- 
ial ZnO films grown on AI2O3 due to the large lattice 
mismatch [t^. Kaiser et al. demonstrated that in high 
fluence Er implanted SiC, the defects act as nucleation 
sites in the formation of Er-atom cluster and NCs [7^ . A 



similar effect can be present in the case of Fe implanted 
ZnO films. 



C. Annealing behavior depends on the initial state 

Note that three kind of samples have been annealed 
at the same temperature of 823 K. One is the ZnO sin- 
gle crystal implanted at 623 K, in which Fe NCs have 
already been formed in the as-implanted states. One is 
the ZnO single crystal implanted at 253 K, in which no 
Fe NCs could be detected in the as-implanted sample. 
The last one is the ZnO thin film implanted at 623 K, in 
which rather small Fe NCs have been formed in the as- 
implanted sample compared to that in the single crystal. 
The annealing behavior in the first case is quite differ- 
ent from the later two cases. One reason is the fact that 
diffusion of Fe NCs is much more difficult than of single 
Fe ions. In the 623 K implanted sample, Fe NCs have 
already been formed, and they are not so mobile during 
823 K annealing. Therefore the size and the amount of 
Fe NCs only slight increase after annealing. However, in 
the sample without Fe NCs or with very small NCs, Fe 
ions are more mobile with annealing, and they aggregate 
into rather larger Fe NCs. Another reason could be the 
same as discussed in the above section, given the fact 
that implantation at 253 K induces more point defects 
than that 623 K (see Figure 2). 



D. Magnetic couphng of dispersed ionic Fe 

Although a part of the implanted Fe ions aggregated 
to metallic NCs, the remaining are in the ionic state. 
Even after 823 K annealing, there is still a considerable 
amount of ionic Fe. By SR-XRD no crystalline Fe-oxides 
could be detected. Therefore, these ionic Fe could be 
diluted inside ZnO matrix. Moreover, Ref. [sO] demon- 
strated that the implantation-induced electrical isolation 
of ZnO is removed after annealing between 773 to 873 K. 
Therefore, the carrier concentration is comparable with 
the virgin sample after annealing at 823 K. However, as 
measured by SQUID down to 5 K and CEMS at room 
temperature, Fe^''" and Fe'^+ are not ferromagnctically 
coupled. 
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In addition to conventional thermal annealing which is 
an equilibrium process, a nonequilibrium annealing tech- 
nique, i.e. flash lamp annealing at a pulse length of 20 
ms, was also used by us [zHl- For an intermediate light 
power, the implantation-induced surface defects could be 
removed without creation of secondary phases within the 
implanted region. However, there is still no detectable 
ferromagnetic coupling between these dispersed Fc ions. 

Moreover, currently the absence of ferromagnetism in 
transition metal doped ZnO is a universal problem. Sev- 
eral groups have shown that transition metal ions, e.g. Fe 
[Zi, Mn [i3[zil, and Co [l3[i3[zil, are substitutional 
inside ZnO. However, no ferromagnetism could be ob- 
served due to the possible reason of the lack of p-type 
conductivity. 

V. SUMMARY AND CONCLUSIONS 

(1) In general, a combination of SR-XRD, ZFC/FC 
magnetization and element specific spectroscopy mea- 
surements is a reliable approach to clarify the observed 
magnetism in DMS materials. 

(2) By correlating the structural and magnetic prop- 
erties of all investigated samples, it is clear that ferro- 
magnetism is only observed when a-Fe (or ZnFe204) 
NCs are present. In as- implanted and 823 K annealed 
samples, dispersed Fe^+ and Fe^+ are the predominant 
charge states. However, they are not ferromagnetically 
coupled. 



(3) a-Fe (bcc) NCs are not crystallographically ori- 
ented inside ZnO matrix. However, fcc-ZnFe2 04 NCs 
formed after annealing at 1073 K are epitaxially embed- 
ded in ZnO. This is due to the crystalline symmetry. 
Hexagonal ZnO crystals are six-fold symmetric, while a- 
Fe is four-fold symmetric. Fcc-ZnFe204 is also six-fold 
symmetric viewed along [111] direction. 

(4) The magnetic properties of these Fe NCs were care- 
fully investigated regarding their memory effect and mag- 
netic anisotropy. A memory effect is observed in the tem- 
perature dependent magnetization measurement, which 
is induced by the different relaxation times originating 
from the different grain sizes of the Fe nanoparticles, and 
consequently different anisotropy energy barriers. The 
in-plane magnetic anisotropy could be due to the shape 
effect. 

(5) The phase separation, i.e. the formation of metal- 
lic Fe, depends on the initial state of the host materi- 
als, namely the carrier and/or the defect concentrations. 
The n-type carriers could facilitate the self-organization 
of metallic Fe NCs. 

(6) The next question is directed to the magnetically 
activation of the diluted ionic Fc in ZnO. The realization 
of p-typc doping for ZnO could be the solution. 
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